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00 ' 

■ We show that a D3/D7 system (at zero quark mass Hmit) at finite isospin chemi- 

o : 

^ , cal potential goes through a superconductor (superfluid) like phase transition. This 
> ' 

O ' is similar to a flavored superfluid phase studied in QCD literature, where mesonic 

^ ' operators condensate. We have studied the frequency dependent conductivity of the 

, condensate and found a delta function pole in the zero frequency limit. This is an ex- 

' ample of superconductivity in a string theory context. Consequently we have found 

Q^' a superfluid/supercurrent type solution and studied the associated phase diagram, 

j-^i The superconducting transition changes from second order to first order at a criti- 

^ ! cal superfluid velocity. We have studied various properties of the superconducting 

O ' system like superfluid density, energy gap, second sound etc. We investigate the 

0^ ■ possibility of the isospin chemical potential modifying the embedding of the flavor 

rn ; 

. branes by checking whether the transverse scalars also condense at low temperature. 

^ \ 

00 ' This however does not seem to be the case. 

o : 

> ' 

X: 

^ ■ I. INTRODUCTION 

The phase structure of QCD (and similar theories) at finite isospin and baryon chemical 
potential is an interesting issue. At high enough isospin density the color singlet mesonic 
flavor degrees of freedoms (e.g. pions) may go through a Bose-Einstein condensation. The 
physical motivation to study such pion superfluid formed at high isospin density is related to 
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the investigation of neutron stars, isospin asymmetric nuclear matter and heavy ion colhsions 
at intermediate energy. Unfortunately this set of problems are difficult to tackle numerically 
due to the complex nature of the action. Various approaches including lattice simulations 
are used to investigate the nature of the QCD phase diagram at finite isospin chemical 



potential and existence of a superfiuid like state is argued 

One way to address various facts about gauge theory is to use the gauge gravity duality 
0] and study a supergravity/tree level string theory to learn about large-N gauge theories. 
Although such examples do not include QCD or even pure YM theory yet, many qualitatively 
similar models has been constructed. At the idealized limit where the ratio of fiavor and 
color degrees of freedom is small, one can introduce probe branes in the gravity background 
to study fiavor physics jsj. In this scenario, the baryon/isospin chemical potential maps 
to chemical potential for various gauge fields living in the brane. The issue of baryon and 
isosgin diemical potential has been addressed in various type of brane systems j^, Q, 111, 



12 



14 



15| . One issue which is relatively less discussed in string theory literature is 



the issue of fiavor superconductivity. There are isospin charged bosonic states with mass of 
0(1) (e.g pions in QCD), which may be thought as strings which has endpoints on different 
fiavor branes. Such a state may naturally condense as we turn on the isospin chemical 
potential. Here we discuss such a scenario in a D3/D7 system in the zero quark mass limit. 
We introduce a couple of branes with no separation between the branes and the branes end 
on a AdS^ black hole. We turn on a chemical potential corresponding to the SU{2) isospin 
gauge field living in the world volume of the branes and study the resulting superconducting 
phase transition and various properties associated with it. 

The plan of the paper is the following. In section [TTl we will set up the probe brane config- 
uration and the equations of motion for the gauge fields. In sections lllll and HVt we establish 
the superconducting phase transition and study the frequency dependent conductivity which 
has a pole (corresponds to DC superconductivity) at zero frequency. The speed of second 
sound is also calculated. Section |V] is devoted to DC supercurrent and the phase diagram 
as a function of the temperature and the velocity of the supercurrent. We will discuss in 
detail the possibility of other relevant adjoint scalar fields also condense at low temperatures 
through a similar mechanism in section IVII In section IVIIl we will conclude and point out 
possible future extensions to the project. 

Note added: When our work is near completion, another paper |l^ came which deals 
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with similar questions. These authors have however kept the full DBI action and the resulting 
details are therefore a httle different. 



II. GENERAL SETUP 

Let us consider AdS^ x in Poincare co-ordinate, 

ds^ = -f{r)dt^ + J^ + + d4 + d4) + L^dfll (1) 

where 

fir) = ^ - ^ (2) 

and L is the radius of the anti-de Sitter space and M is the the mass of the black hole. In 
this note we will adopt the convention that M — L — 1. The temperature of the blackhole 
(and also the boundary field theory) is given by 

r=i. (3) 

TT 

It is more convenient to analyze the system by making a coordinate transformation z — 1/r. 
The metric becomes: 

ds^ = -f{z)dt^ + -^j— + -^{dxl + dxl + dxl) + dnl (4) 

with 

/(-) = 4 - (5) 

The horizon is now at 2; = 1, while the conformal boundary lives at 2; = 0. In this back 
ground we will introduce two (or possibly more) co-incident D7-branes. For simplicity we 
will consider the zero quark mass embedding, where the brane is filling the whole AdS'^ and 
wrapping the maximal of . In this limit the effective induced metric on the brane will 
just be the AdS^ x S"^ metric. 

dz^ 1 

^«Lne = + + -^{K + + <i^\) + <i^\ (6) 

The effective action for the brane fields is the Born-Infeld action, 

Sj^jij = -T7 J d^xVG + 27ia'F (7) 
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In order to consider the system at finite isospin chemical potential, we will put in a pair of D7 
as probe branes. In this case, F is an f/ (2) field strength on the world volume of the probes. 
We will not focus on baryonic U{1)b and only investigate terms containing SU{2) isospin 
gauge fields. The string states which has their endpoints at different D7-branes are charged 
under the isospin SU{2). The exact form of the non-Abelian DBI action is unknown [17]. To 
proceed further we will expand the action to leading order in Xyma keeping only Yang Mills 
terms. Such simplification has also been employed in studying other aspects of holographic 



QCD such as the meson spectra and baryon masses [18|, ll9|]. Such an approximation will be 
more accurate in the limit where the non-Abelian field strengths are small, which may not 
be the exact situation in our case. However we believe that essential physics is captured. 
The effective action now takes the form 

Sj,Bj = -Ty^^^^ J (fxV^TrF^ (xN.J d'xV^TrF^ (8) 

where we have scaled out 7+1 dimensional Yang Mills coupling gj and 

F» = dA^ + e^^^A^A^ (9) 

The setup is very similar to that of 



20|, where the non-Abelian gauge field is shown to 
condense at low temperature in an AdS black hole background. Due to the non-Abelian 
nature of the SU{2) symmetry, the and components of the gauge field are charged 
under the component. Hence turning on a chemical potential for the component of 
the gauge field may lead to a condensation of the other two. This mechanism is very similar 
to condensation of U{1) charged scalar discussed in [2]]], where a charged black hole may 
lower its free energy by distributing its charge in a scalar condensate. Turning on a chemical 
potential for the component of the gauge field breaks SU{2) to f/(l)/. A condensation 
of or component of the gauge field further breaks the U{l)j symmetry. It should be 
mentioned that unlike U{1)b where all charged states are baryonic and has a mass of order 
0{N), this U{l)i theory has charged states (mesons) with 0(1) mass and their condensation 
can naturally be studied in terms of the probe brane picture. 
We start with the ansatz 

A = AtT^dt + B^T^dxi (10) 

We will assume spatial homogeneity in the field theory directions and our fields will only 
have dependence on the radial coordinate. The equations of motion for the fields in this 
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coordinate system are: 



4' - y - 2||a = (11) 



B:+iij + z]B'^ + ^iA't^) =0 (12) 



To require regularity at the horizon we will have to set = at 2; = 1. Since we have a 
set of coupled equations, this will in turn give the following constraints at the horizon 

B'. = 
At = 

at 2; = 1. Examining the behavior of the fields near the boundary, we find 

At ~ iJ, — pz^ + ... 

Using gauge/gravity duality, /x, p are mapped to the isospin chemical potential and charge 
density in the dual field theory, respectively. is mapped to the expectation value of a 
meson operator which condenses at low temperatures. We will set the non-normalizable 
mode Mr to zero. In the following we first establish that the mesonic condensate forms 
below a critical temperature. We compute the time dependent conductivity by turning on 
a spatial component for the isospin current as fluctuations. 

A^ = X{z)e'''^T^dxz (13) 

The equation of motion for Aj. is 



We will choose ingoing boundary condition at the horizon A^ oc {z — 1) '^1^. Asymptotically, 

Xr^S^^ J^Z^ + ... (15) 

Jx corresponds to the isospin current, while gives the dual current source (superfiuid 
velocity). The conductivity is given by 

J X 



Re 



iujSx 



(16) 
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X can be normalized to one at the horizon. As we will see, the conductivity has a pole at 
uj = 0. This suggests that there is a DC supercurrent solution. To find such a solution we 



solve the following set of coupled time independent equations 22|, |23|]. Note here the effects 
of Ax = X{z) on the other components of the gauge fields are taken into account, so this 
is not a fiuctuation around the condensate formed by Bx. In this case the field A^ has the 

2 



form A^ = X(z)T^dx'i. 



A' B 

A'; -^At = 

z z^f 



Uj + -]b'^ + -^ A^,^-BI-X^bA = (17) 



X" + (i + L)x' - Bl^=0 
The consistency condition At = at the horizon also imposes the following constraint on 



Ax at the horizon 



X' = - Bl^ (18) 



As we will see shortly, the system Eq. (1171) reveals a rich phase structure as the boundary 
value of Ax is tuned. 

The convenient physical parameters for us are (^; ^^^^ ), or(^, ^, "^^^~)- We 

will use -^^^-is an order parameter and plot the the phase diagrams as a function of (^, ^). 
In practice, we choose to keep the temperature fixed and vary /i in this paper. The com- 
ponents of gauge fields on the three sphere may also condense through similar mechanism. 
We will examine these cases in the appendix. It is also an important question to check 
whether the isospin chemical potential would modify the embedding of the fiavor branes. 
We will leave the detail of these investigations in a future work. Here we will begin look- 
ing into these questions by checking whether the transverse scalars also condense at low 
temperature (section IVip . 

III. PHASE DIAGRAM 

At high temperature (equivalently small /x) there is only one branch of solution to the 
set of equation Eq. (fT2|) given by. 

At = Ml-^') (19) 
Bx = 
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FIG. 1: Plot of the zero mode at /i = 4 



This should be interpreted as isospin-charged black hole, where gauge fields are confined in 
the D7 brane. The dual gauge theory interpretation is a deconfined plasma with non-zero 
isospin charge. As we increase /i the effective mass of Br^ in Eq. (fT2|) decreases and 
develops a zero mode at = = 4. The existence of this zero mode can be analytically 
demonstrated. Putting the value of At from Eq. (fT9|) in the second equation of Eq. (fT2|) and 
we get (this small fluctuations analysis does not depend on the any possible cubic terms and 
hence true for other scalar fleld ansatzs considered later) , 



The plot of the zero mode is shown in Fig[T], 

Any further increment of /i leads to a condensation of B^. Hence for < 0.25 there 
is a new branch of solution with non-zero value of B^. Such a solution can be numerically 
constructed. The associated transition seems to be a second order transition from our 
numerics. Here we show the plot of condensate (the plot of free energy is provided in section 



We flnd that the condensate levels off at low temperature to W^/fJ^^ ~ 0.26. In terms of 
more familiar variables critical temperature (Tc), the condensate strength can be expressed 




(20) 



The above equation has an analytic solution for = 4, given by, 




(21) 



EH) with 



as ^ 0.264^7T^T^ ^ 515.94T3 or VT^ ^ S.OIT^. 



A. Speed of Second Sound 



The boundary fleld theory which is dual to the D3/D7 system in AdS^ x behaves like 
a superfluid below the critical temperature. Superfluids are known to exhibit modes known 
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FIG. 2: Plot of the condensate with 1/// 




0.15 



0.20 



0.25 



FIG. 3: Speed of second sound as a function of 1/// 



as second sound which are basically temperature waves propagating through the fluid. For 
a hydrodynamic discussion see [23|, where this was computed at zero superfluid velocity for 
the Abelian Higgs model on AdS^. We also compute the speed of second sound in our case. 
The superfluid velocity now corresponds to Sx/ fJ., which we set to zero for this computation. 



The main relation is (see Eq. (18) of 23|): 



Ps 



vl = (22) 

where ps is the density of the superfluid component and P is the pressure. The pressure can 
be expressed in terms of the total fluid density p by using the equation of state of a perfect 
fluid P = p/{d — 1) where d is the dimension in the fluid (boundary) theory^. Using this it 



^ The perfect fluid approximation is valid here since we are not considering any backreaction due to the 
metric. Hence there is no viscosity correction which originates from fluctuations of the metric. 
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Im[-] 




(a)Plot of real part of - 




(b)Plot of imaginary part of - 



FIG. 4: Plot of the real part of ^ with /i = 7.57/Zc, 2.52/Xc, 1.71/Xc, 1.37;Uc, 1.13;Uc (gradually from 
red to green curves). The blue curve is for the exact frequency response at ^ = = 4. 

is fairly simple to compute f|. We present the result in Fig [31 where we plot f| as a function 
of l//i. At high values of /x f| approaches a limiting value f| = 0.013. 



IV. FREQUENCY RESPONSE 



Here we will study the frequency dependent conductivity of the spatial component of the 
isospin current Eq. (fT3!) . Without any condensation of the frequency response can be 



exactly solved 



24| , and is expressed in terms of digamma functions. With the condensation 



of field, we can still numerically calculate conductivity using Eq. fITB]) . We solve the 
equations Eq. f[T^ to get the condensate and on the fixed background given by the solution, 
we solve Eq. f|T^ . The corresponding a is shown for various values of the parameter n 
in Fig m We find that Im[a{uj)] ~ ^ as — 0. Where is the superfluid density. 
The plot of the superfluid density is presented in Fig [51 Near fi = fi^, Ug vanishes linearly 
and is proportional to /i — fic- Fitting a linear function near the critical point we get 
Us oc O.lfiKfi — He). A pole at = implies a i?e[(j(0)] ~ TinsSiuj) + terms regular in uj. 
This delta function singularity of the real part of sigma is not captured in the numerics 
directly. However this corresponds to the superconductivity /superfluidity and consequently 
we can find a supercurrent/superfiuid type of solution (see Sec [V|). Unlike 1^ we do not 



10 




FIG. 5: Plot of superfluid density with 



get any low temperature resonances in the conductivity. Our result is more similar to the 

2 



zero mass Abelian-Higgs system presented in \24i\ . 



As J Re[cr]da is a temperature invariant quantity, the delta function at = is com- 
pensated by a dip of Re[a] for lower value of frequencies. Also the curves for various 
temperatures approach the no-condensate curve for higher frequencies. The dip becomes 
more prominent as we lower the temperature (increase fi). It is clear from the diagram that 
at low temperature (for large /i) i?e[cr] 0. In fact it is expected that at low temperature 
Re[a] ~ exp(— ■^), where Ag is the energy gap of the system. Also looking at the zero 
temperature limit of the real part of conductivity we see that i?e[ti;] = for < Ap. Ap is 
similar to the energy of a "Cooper pair". The ratio n-g = ^ gives important information 
about the nature of the condensate. From our numerics we calculate, 

Ug ^ 1.2 (23) 



V. EFFECT OF STATIONARY ISOSPIN CURRENT 

In this section we investigate the effects of turning on a finite time-independent isospin 
X field (recall ^^(r) = At{r)T^dt + Bx{r)T^dxi + X{r)T^dx3). We tune the boundary value 
of X, Sx/fJ' at a fixed isospin chemical potential fi to investigate the phase structure. At 
high enough /i and in absence of X the gauge component condenses. As we can see from 
Eq. (fT7|l . the effect of X is to increase the effective mass of the Bx field. This will cause the 

^ It seems that in an Abelian-Higgs system in AdS^ resonances occur near conformal mass. [2^ 
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FIG. 6: Plot of Wx/fJ-^ as a function of Sx/ fJ- at different values of n- The values of /i range from 
1.24/ic to 1.97//C from left to right. 

Bx condensate to weaken with increasing Sx/ /i- Indeed, this happens, as we can see from Fig 
ini Here we plot the condensate strength as a function of the isospin current source Sx/^ for 
different chemical potentials /i. fi increases from left to right. For strong enough Sx/ ^ (above 
a critical value) there is a phase transition to the normal (non-superfluid) state. The order 
of this phase transition seems to be /i dependent. For high /i (compared to /ic) the phase 
transition is first order, i.e. the system discontinuously jumps to the normal state above the 
critical current velocity Sx/ /i- For /i close to fic the transition becomes second order. The 
order of the transition changes near fi^p = lAfi^- Note that in each case the condensate 
approaches a limiting value at low values of Sx/fJ', and this limiting value decreases with 
decreasing /i. 

In order to properly see the transition from first to second order, we also plot the difference 
in free energies of the normal and superconducting branches. In Fig [7] the left hand plot is 
for yU = 1.97yUc while the right hand plot is for fi = 1.24/ic. We see the typical swallow tail 
shape indicating the abrupt change in dominance from the normal to the superconducting 



branch at low values of Sx/l^- We can understand this curve Fig 7(a) in the following 
manner. As Sx/ /i is lowered from above the critical value Sx^c/ f^, at first there is only a 
normal branch (I). At some value Sx,N/f^ = 0.571 two new branches are nucleated: one 
of these is stable (II), while the other one is unstable (III). The stable branch starts out 
with a higher free energy than the normal branch, but as Sx/fJ^ is lowered further these two 
branches intersect at Sx,c/ 1^^ = 0.568, where the branch (II) has the same free energy as 
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FIG. 7: Plot of the difference in free energies of the normal and superconducting branches as a 
function of S^^Z/i 

branch (I). This is the first order phase transition point below which the system jumps to 
the superconducting branch (II) which now has lower free energy. At lower values of the 



chemical potential shown in Fig 7(b) the transition is continuous between branches (I) 
and (II) at S^J^i = 0.038. 



VI. WHAT CONDENSES AND WHAT DOESN'T 

There are various type of bosonic fields living in the branes which are charged under 
and may condense as we turn on a chemical potential for A^. These states come from strings 
ending at different branes. This include the example of gauge fields proportional to and 
r^. Together with A^, the trio corresponds to a vector iso-vector operator (flavor current) ^ 
in the boundary theory. We have considered a ansatz like Eq. ffTOl) . but more generally one 
may consider A = AfT^dt + B^r^dxi + CxT'^dx2, then action contain a repulsive interaction 
term (oc -B^C^) between the two, so one may expect that both do not condenses together. 
And a ansatz of the form B^lr^dxi + T'^dx2) t20i] will generally have a higher free energy 
than the ansatz we have chosen. 

Another possibility is the transverse scalar fields in the brane. Brane fields are invariant 



As niq = to begin with in our case, there is no stable meson. We only have quasi normal modes 
corresponding to various brane fields 26[ . The dual interpretation of such modes are mesons decaying in 
gauge theory plasma. 
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FIG. 8: Plot of action for phases Bx ^ (upper curve) and </> / (lower curve). 



under only a 50(4) x 5*0(2) subgroup of the full R-symmetry group 50(6) of 5^. Two 
transverse scalar fields correspond to one 50 (2)-charged iso- vector scalar operator in the 
boundary theory. The Born-Infeld type effective action for such scalars is given in the 
appendix. We choose a general ansatz suitable for our case, 

$1 = 0ir\ $2 = (f)2T\ (24) 

The EOM's are, 

<P: + + ^W. + ^ (f - = 0' (25) 

where i,j = 1,2 and i ^ j From the discussion of the previous paragraph it is clear that 
due to the repulsive interaction term, a preferred configuration will have either 0j's turned 
off. One may try to find when such a field becomes unstable in a fixed At background. 
From our numerics we find out that for the solution Eq. fl20l) such a thing happens at 
/i^ ~ 6.57 = 1.64/ic and greater than our /i^ = 4 value for field. Hence when we gradually 
decrease the temperature of our system, condenses before any scalar degree of freedom. 
One may further ask wheather the resulting superconducting phase with a B^. condensate 
has a instability for fluctuations. Our numerics answer this question negatively. It seems 
that some type of "blocking" mechanism stops further condensation of more fields. However 
that does not rule out the possibility of a first order transition between B^ condensed phase 
and (f) condensed phase. We plot the free energy of both the phases to investigate a possible 
first order transition Fig [HI It seems that the B^ ^ Q phase always dominates. However 
our numerics is not good for the parameter range /i > 10/ic, also we did not exhaustively 
searched for the possibility of various mixed phases. We also have to address the issue of 
curvature of 5^ in more detail. The case for gauge fields with 5^ indices is also similar to the 
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scalar fields (see appendix), at least for small fluctuations. Hence the blocking mechanism 
discussed before will work for them too and these fields will not also condensate. However 
we have not investigated the possible phases and possibility of a first order transition in 
great detail for these fields. 



VII. CONCLUSIONS 



In this paper we explored a model which embeds the mechanism of holographic super- 
conductivity/superfluidity in a string theory setup. We have studied a couple of probe D7 
branes in AdS^ x back ground. We turn on the finite isospin chemical potential (i.e. 
potential for some of the world volume gauge fields) and have found the existence of a fla- 
vored superconducting state at high enough chemical potential. We studied the frequency 
dependent conductivity and found a delta function pole in the zero frequency limit. This 
is superconductivity. Consequently we have found a superfluid/supercurrent type solution 
and studied associated phase diagram. The superconducting transition changes from sec- 
ond order to first order at a critical a superfluid velocity. The holographic dual of such 
string theory system is large N^, A/" = 4 supersymmetric gauge theory with Nf -C A^^^. In 
a dual gauge theory such a superconducting state is characterized by mesonic condensates. 
We have studied various properties of the superconducting system hke energy gap, second 
sound etc. 

Here we have discuss a possibility of a first order transition between various possible 
condensate. It is also an important question to check whether the isospin chemical potential 
would modify the embedding of the flavor branes and whether the transverse scalars may 
condensate. In our case it does not, for some other setup it might. It would be interesting 
to address such questions in more detail. 

In QCD, the mesonic condensates (pion superfluid) have been argued to exist at finite 
isospin chemical potential. A natural extension of our work will be to study more realistic 
holographic models like 1^, Q or in the context of AdS/QCD like theories. It would 
be interesting to study such phenomenon where the fundamental degrees of freedoms are 
fermionic. Here we also focused on the zero quark mass case, which may be thought of 
as a high temperature limit. It would be natural to extend our work to the case of finite 



quark mass. Another interesting question to investigate is the flavor backreaction 



271, 



28|, 
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29|. The mesonic operators which condensate in our case are color singlets and do not 



lead to color superconductivity. It would be interesting to study some models with color 
superconductivity. 
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Appendix 

Consider D7-branes embedding in AdS^ x background, 
dz^ 1 

ds'^ = -f{z)dt^ + ^ + —{dxl + dxl + dxl) + dO'^ + sin^ 6 {d^p^ + sin^ Lp dnl) , (26) 
the induced metric is 

dz"^ 1 

ds"^ = -f{z)df + ^Y[z) ^ i2 ('^^i + + dxl) + sin^ ^ sin^ d^l. (27) 



According to [XT'], the leading order DEI action of Dp-brane is 

= -^^^^^^ j d^+^i y^G^^Tr [F,,F^' + 2V,^{D'^^' + ^^]] , (28) 

where scalars = X*/ (27ra') are the transverse coordinates, and covariant derivative is 
defined as Va^i = da^i + K, <l>i]. 

For D7-branes, one has two scalars i = 6,(f. 

If we turn on the gauge field with ansatz, A = AtT^dt, where = T"'/2i, and thus has 
the commutation relations [T",r^] = e°''"^T'^. The scalars take the general form, = ^^r". 
The non-zero components of 'Da^i are 
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The effective action can be written as 



T-{2^(/f /• 8 siii^ ^ siii'V 



Tr {FabF'^') 



1=1 



- sin^ e 

2 



1,2,3 



V\2 



(30) 



wliere ga are 10-dimensional metric witfi i — 9 or (p. 

Introduce = {^\ ± ^$2)/2, we find tfie equations of motion for ($3, $^), 



9ii 



gu 



fiz) 



z 



z^ z 



1 

is 



$V + ^'-(-(H)^ + (^3)^) - *3 ^-^3 



Z^ z 



9 it a2 



where i,j G {6, ip) and i 7^ j. 

Consider fluctuations, one can set sin^ = 1, and thus ga = 1 in above equations. If we 
restrain our discussion in the ansatz 



the EOM's can be simphfied, 



0, 



(32) 



(33) 



where j = 1, 2 and i ^ j. 

For the simplest case, with $ = (f){z){T^dx^ + r'^dx'^), the equation of motion is 



(34) 



If the gauge fields on the 3-sphcrc instead of the position scalars of D7-branes are turned 
on, one would find similar condensation. In an ansatz 



A Atiz^dt + Ae{zyde, 



(35) 
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the equations of motion are 



Ae+{-Z + ^]A'o + -A-Ae = ^- (36) 



Or in another ansatz, 

A ~ AtT^dt + Ae{T^de + sin Odip), (37) 

where {9, (p) are angle coordinates on 3-sphere of D7-brane worldvolume AdS^ x S^. The 
EOM's turn to be 



[1] K. SplittorfF, D. T. Son, and M. A. Stephanov, Phys. Rev. D64, 016003 (2001), hep- 
ph/0012274. 

[2] D. T. Son and M. A. Stephanov, Phys. Atom. Nucl. 64, 834 (2001), hcp-ph/0011365. 

[3] D. T. Son and M. A. Stephanov, Phys. Rev. Lett. 86, 592 (2001), hep-ph/0005225. 

[4] D. Toublan and J. B. Kogut, Phys. Lett. B564, 212 (2003), hep-ph/0301183. 

[5] J. B. Kogut and D. K. Sinclair, Phys. Rev. D66, 034505 (2002), hep-lat/0202028. 

[6] X. Hao and P. Zhuang, Phys. Lett. B652, 275 (2007), hep-ph/0612320. 

[7] J. M. Maldacena, Adv. Thcor. Math. Phys. 2, 231 (1998), hep-th/9711200. 

[8] A. Karch and E. Katz, JHEP 06, 043 (2002), hep-th/0205236. 

[9] M. Rozah, H.-H. Shieh, M. Van Raamsdonk, and J. Wu, JHEP 01, 053 (2008), 0708.1322. 

[10] D. Mateos, S. Matsuura, R. C. Myers, and R. M. Thomson, JHEP 11, 085 (2007), 0709.1225. 

[11] J. Erdmenger, M. Kaminski, P. Kerner, and F. Rust (2008), 0807.2663. 

[12] S. Kobayashi, D. Matcos, S. Matsuura, R. C. Myers, and R. M. Thomson, JHEP 02, 016 

(2007), hep-th/0611099. 

[13] J. Erdmenger, M. Kaminski, and F. Rust, Phys. Rev. D77, 046005 (2008), 0710.0334. 

[14] O. Aharony, K. Peeters, J. Sonnenschein, and M. Zamaklar, JHEP 02, 071 (2008), 0709.3948. 

[15] K. Ghoroku, M. Ishihara, and A. Nakamura, Phys. Rev. D76, 124006 (2007), 0708.3706. 



18 



[16] M. Ammon, J. Erdmenger, M. Kaminski, and P. Kerner (2008), 0810.2316. 

[17] C. V. Johnson (2000), hep-th/0007170. 

[18] T. Sakai and S. Sugimoto, Prog. Theor. Phys. 114, 1083 (2005), hep-th/0507073. 

[19] H. Hata, T. Sakai, S. Sugimoto, and S. Yamato (2007), hep-th/0701280. 

[20] S. S. Gubser (2008), 0803.3483. 

[21] S. A. Hartnoll, C. P. Herzog, and G. T. Horowitz, Phys. Rev. Lett. 101, 031601 (2008), 
0803.3295. 

[22] P. Basu, A. Mukherjee, and H.-H. Shieh (2008), 0809.4494. 

[23] C. P. Herzog, P. K. Kovtun, and D. T. Son (2008), 0809.4870. 

[24] G. T. Horowitz and M. M. Roberts (2008), 0810.1077. 

[25] Current authors: Work in progrss. 

[26] P. K. Kovtun and A. O. Starinets, Phys. Rev. D72, 086009 (2005), hep-th/0506184. 

[27] R. Cascro, C. Nunez, and A. Paredes, Phys. Rev. D73, 086005 (2006), hep-th/0602027. 

[28] S.-J. Sin, JHEP 10, 078 (2007), 0707.2719. 

[29] P. Basu and A. Mukherjee, Phys. Rev. D78, 045012 (2008), 0803.1880. 



